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Abstract. One-dimensional scattering by a Coulomb potential V{x) = is 
studied for both repulsive (c > 0) and attractive (c < 0) cases. Two methods 
of regularizing the singularity at x = are used, yielding the same conclusion, 
namely, that the transmission vanishes. For an attractive potential (c < 0), 
two groups of bound states are found. The first one consists of regular (Rydberg) 
bound states, respecting standard orthogonality relations. The second set consists 
of anomalous bound states (in a sense to be clarified), which always relax as 
coherent states. 
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1. Introduction 

One-dimensional quantum Hamiltonians are very useful for modeling simple quantum 
systems. Beside their ubiquitous importance in the study of transmission and 
tunneling experiments, numerous quantum systems in higher dimensions can be 
reduced to one-dimensional ones, due to symmetry (for instance radial wave functions 
in a central potential) or specific physical properties (Josephson junctions or edge 
states in the quantmit Hall effect arc just two examples). 

The aim of the present work is to examine one-dimensional scattering by a three- 
dimensional coulomb potential V{x) = j^^^j^, starting from the Schrodinger equation 

with Hamiltonian H = ^ + V, for an eigenstate ip{x), with a; G M* = M \ {0}, 

-§(^) + ^^(a^)=e^(^), (1) 

with A = ^""^^2 and e = where E is the energy. A > corresponds to the repulsive 
potential, A < to the attractive one; The boundary conditions will be specified later 
on. This is referred to as the one- dimensional Coulomb potential problem. Although 
it has recently been studied[T], we find it useful to analyze it using somewhat different 
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approach. As it turns out, there are some subtleties involved, which might affect some 
of the conclusions reached in Ref. [T] . 

One of the main advantages encountered in the quantum Coulomb problem is 
that the exact wave functions are computable. In three dimensions, it has been shown 
eighty years ago [2] that the asymptotic behavior of the wave functions is somewhat 
distinct from that of plane waves. This property has been shown to be valid also in 
one dimension [3j. 

It proves useful to follow, first, the standard reduction of the Coulomb problem 
in three dimensions into a radial one-dimensional equation, and to point out the 
differences between this equation and Eq. ([T]). Starting from the three-dimensional 
Schrodinger equation, carrying out partial wave expansion 'l'(r) = Si^o(2' + ^) 
ip I (r) Pi {cos 9), and writing the radial wave function as V'i(^) = one obtains 

the radial Schrodinger equation for 4>i{r), with < r < oo, 

l{l + l) A" 



dr^ 



Mr)^eMr)- (2) 



For I = {s wave scattering), Eq. ^ has the same form as Eq. ([I]). The two basic 
solutions of Eq. ([2]) are the regular one, satisfying 0z(O) = 0, and the singular one, 
satisfying (^i(O) ^ 0. The singular solution should be discarded: if not, for / > 0, the 

probability of finding the particle in a sphere of radius R, Pi{R) = pi{r)2Trr'^dr 
becomes infinite for any R ; for I = 0, the situation is more subtle, Po{R) remains 
finite, but the radial current JoiR) = /q jo{r)2Trr^dr > becomes nonzero, which is 
impossible for an s state [31[S]. 

A couple of difficulties arise when Eq. ^ is considered as compared with Eq. ^ : 

(i) The solutions of Eq. ^ are required on R* , and not only on its positive part R*^ . 
Note that H is invariant under space inversion. 

(ii) The arguments used in the three-dimensional case to discard singular solutions 
of Eq. (21) are not valid [5] for the original problem specified by Eq. ([1]), and 
the imposition of scattering boundary conditions requires their inclusion as well. 
The standard techniques used for matching the wave function at .t = require 
either the calculation of ip'ie) or of /_ V{x)dx and both quantities diverge 
logarithmically when e — > 0. One must then cope with ultraviolet divergences, 
which need to be regularized. 

These difficulties lead us to the connection problem, which can be defined as follows: 
Let us decompose Eq .([1]) into two equivalent coupled equations, one defined on W*^ 
with V{x) = ^, the general solutions of which read 

tp+ix) = Af{kx) + Bg{kx) , (3a) 

and the second defined on Ml with V{x) = — ^, the general solutions of which read 

?/>_ (x) = af{kx) + bg{kx) . (3b) 

Here, /(x>0) and f{x<0) are regular solutions, while g{x>0) and g{x<0) are singular 
solutions, defined on the appropriate domains; the relations between f,g and f,g will 
be clarified later on. The connection problem consists in the calculation of the 2x2 
matrix expressing (A, B) in terms of (a, b). Since the derivative of the singular solution 
diverges at a; = 0, it is impossible to match both ■0 and ip' at x = 0. It is also not 
possible to use the method [71 [8] employed in a problem of scattering by a potential 
V{x) = \6{x) since the latter potential is integrable at a; = 0, J^^V{x)dx = A, whereas 
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the Coulomb potential is not. Apparently, the connection problem cannot be solved 
in terms of simple linear relations, and one needs to consider bilinear constraints 
(an example of such a constraint is the current conservation J(0~) = J{0~^) around 



Our first task is to properly formulate and solve the scattering problem, 
corresponding to e > 0. To carry it out, we use two independent regularization 
methods. One is based on bilinear constraints, which can be formulated in such a way 
that ultraviolet divergences are canceled. The other method consists in calculating 
the exact transmission for a truncated Coulomb potential Vg, with T4(x) = for 
|a;| < e, Ve{x) = \/\x\ for |a;| > e and letting e ^ 0. With both methods, we arrive 
at the conclusion that the transmission coefficient vanishes, T = 0. The potential 
is perfectly reflective. Moreover, this property of total reflection also holds for the 
attractive potential (A < 0), whereas classically the reflection vanishes; it is a novel 
manifestation of perfect quantum reflection from an attractive potential. It is distinct 
from the standard example of quantum reflection from an infinite attractive square 
well: in the latter case, the divergence of ^V{x)dx is faster than logarithmic, and the 
corresponding spectrum is not bounded from below. 

Our second goal is to calculate bound state energies and wave functions for an 
attractive potential (A < 0) (the one-dimensional "hydrogen atom" problem). The 
ensuing discrete part of the spectrum (e < 0) appears to be rather intriguing, as it 
is composed of two interlacing spectra. The first one (reported also in Ref. [U |9j) 
is the usual Rydberg spectrum, with energies = — with n = 1,2,... The 
corresponding wave functions are the regular solutions of the differential Eq. ([T]). The 
energies of the second part of the spectrum are shifted from the first ones through 
n — > n + 1/2, that is. En = — f^^+Xy ' '^i^h n = 0, 1, . . . The corresponding wave 
functions will be refereed to as anomalous states, and are constructed in terms of the 
singular solutions of Eq. ^ . These solutions are square integrable but not orthogonal. 
A proper incorporation of such states might require further insight into the basic 
principles of quantum mechanics. 

We organize the rest of the paper as follows: In section [21 we will first study the 
scattering problem, then explain, in section [3l the two regularization methods used to 
solve the connection problem. The bound state problem will be analyzed in section [H 
where regular and anomalous states are introduced. Finally, a short discussion of 
our results is carried in section [S] Calculations requiring technical manipulations are 
collected in the appendices. 

2. The scattering problem 

2.1. Scattering states 

2.1.1. Basic solutions For the scattering problem, we have e > in Eq. ([T|). It is 
convenient to recast Eq. (U) so that all quantities are dimensionless. Let k ~ -y/e, 
u = kx, Tj = \/{2k) ~ h \^JE ^^"^ ~ ^(f)- Then the equation for tp is 



X = 0). 




(4) 
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Figure 1. (full line) and (dashed line) for 77 = 1/5. 



with regular and singular solutions frj{u) and gri{u). Eq. ([¥]) is equivalent to the 
following couple of equations : 

- —^{u) + 2-ip{u) ^ ip{u) for u > ; (5a) 

72 

— r^H-2-(/9(u) ^Lp(u) foru<0. (5b) 

The solutions of Eq. ((5a|) arc known as Coulomb s wave functions [2l[T0] with L = 0. 
We will write ^l,(u) the regular solution and Gjj{u) the singular (logarithmic) one: 

F,j (u) = C^u c~'" Af (1 - ir;, 2, 2iu) ; (6a) 



^^ »e ■"!-( ir;) ^^^ _ .^^ 2^ 2iu) - i(-l + vrry + 2i^)F^(?/)/C2 , (6b) 



where 



' *^/^h^^^ ^""^ = r7Q(r(l - ir;)) 



In these equations, AI is the regular confluent hypergeometric function, also written 
as 1^1, and U is the logarithmic (also called irregular) confluent hypergeometric 
function [TT]. Both and Gj, are real. Thus, the solutions of Eq. ^ for u > 
are /^(u) = F^(w) and g^(u) = G^(ri), Vr]. 

Consider now the domain u < 0. In principle, finding the solutions of Eq. (|5bp 
can be achieved by direct continuation of Fri{u) and G^(u). Practically, this requires 
some care, especially for G^. Fj^ can be continued analytically since it is regular at 
u = 0, while for G,,(m) one has to avoid the divergence of G^ at m = 0. Since (|5a|) is 
valid for any sign of r], we simply need to change rj —77 in the previous expressions, 
to get the solutions of (f5b|) . thus we get /?,(«) = F_^(m) and g,^{u) = G-^{u) Vu < 
and V77. It should be pointed out that, in the imaginary part of (|6bp . the factor before 
does not follow the ry ^ — ry transformation [12] . The right expression is (note that 
G_,, = c^'^C^), Vu < : 

g^{u) = -2r; ^ " U{1 + ir;, 2, 2iu) - i(-l + Trr; + 2t„)F_„(w)/G2 ^ . (6c) 

One should also note that relations (14.1.14) to (14.1.20) of [TU] extend for p < as 
soon as one replaces log(2p) by log(— 2p) in (14.1.14). 
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Basic solutions fr]{u) and g,f{u) arc defined on M* and shown on Fig. [T] These 
solutions are constructed so that Eqs. (|5al5bp are satisfied for both u > and u < 0, 
yet the matching condition at m = is not addressed yet. This will be carried out 
when we solve the connection problem. 



2.1.2. The general solution Having defined the basic solutions^ we can now form the 
general solution as a linear combination of fr/iu) and on each side of m = 0. We 
use expressions ((3a|) for u > and (j3bp for u < 0. Now, the relation between / and 
/ and that between g and g are well established, so that bar" can be omitted. With 
these notations, the general solution writes 

1 Af„{u) + Bg^{u) for u > ; 

V[U,T]) = I (7) 

I ajn\u) + ogrj(u) for u < v . 

The linearity of Schrodinger equation implies that the connection problem 
eventually reduces in finding the 2x2 matrix D, which obeys 

s)=d{1) ^ith det{D)^0. (8) 

2.1.3. Transfer matrix It should be stressed that D is not the transfer matrix T 
because T transforms incoming and outgoing (distorted) plane waves at u — > — oo to 
those at u — > oo. In order to identify these asymptotic waves, we need first to examine 
the asymptotic behavior of the function ip{u,r]) when u — !■ ±oo. 

The asymptotic behaviors of F^i^u) and G,j(w), for u +00, have been established 
a long time ago in Ref. [2]: 

Fniu) = (1 + ^ + ^4^+-) - + - ^^^8^^+-^ " ^"^""^^ 

sin(u — B?j(w)) ; (9a) 

cos(u — Qjj{u)) ; (9b) 



u — voo 

with 

e„{u) ^ 7ylog(2w) - arg[r(l + i?/)] . (10) 

Derivation of the asymptotic behaviours of Fri{u) and G,,(u), for u — > — oo, is 
more subtle. Their determination (6c) and (6d) of Ref. [1] is to be reconsidered [13]. 
In Appendix lAl we find 

77 577^ — 77^ 
2^ ^ 8m2 

'2m 8^2 

e-''''sin(u- e^(M)) ; (11a) 

77 577^ — 77^ 
2^ ^ 8u2 
,77^ 277 — 477'^ 



F,{u) = e--^{\ + J- + ^^+..)s\n{u-Q,,{u)) 

+ e--^i^-^^l^+..)cos{u-e,{u)) 



G,{u) = e-''(l + + ^^+..) cos(u ~ e„(u)) 



e-(^-^^+..)sin(.-e,(.)) 



e''''cos(u- e^(M)) . (lib) 
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Thus, the asymptotic form of the solution (p{u, ry), is 
(p{u,ii)'~~^ Asm{u — 6,,(tt)) + Bcos{u — 6,,(m)) 

iA 



(p{u,'q) 



U — > — QO 



(e^(u)-ti) . 

2 ' 2 " 

ae'"' sin(u + e,,(-u)) + 6e~'"' cos(m + e,,(u)) 



6 c-'''' - iae" 



i(6,(u)+u) 



(12a) 



(12b) 



2 2 
The transfer matrix T relates the coefficients of the distorted plane waves at u oo 
with those at u ^ — oo: 



B-iA 
B + iA 



T 



be-'''' - 



(13) 



_6c-^'' + iae'''' 

Solution of the scattering problem is equivalent to elucidation of the transfer matrix. 



2.2. Scattering 

2.2.1. Transmission and reflection amplitudes Alternatively, we define transmission 
t and reflection r amplitudes in terms of a wave (fa propagating from — oo (a = L), 
or from oo (a = R). Explicitly, 



gi(ti+e^(M)) g-i(ti+e^(ti)) 
t-L e''("~®''(")) ; 



and 



-i{u-Qr,{n)) 



me 



i(u-e,,(u)) . 



<y5R(u,?7)< 



u — > — oc 



Time reversal invariance implies t^ ^ ti^ = t and reflection symmetry H{—x) ~ 
H{x) implies tr = tl = r (to demonstrate it properly, one must note that, if (p{u, rj) is 
a solution, tp{^u,r]) is another solution, a priori independent of the first one). Some 
useful relations expressing A, B, a, b in terms of t, r are given in Appendix iBl 

The corresponding transmission and reflection coefficients are 



T^\t\ 



R 



(14) 



and fulfill i? + r = 1 (see Eq. ([39a|). For t ^ 0, it is instructive to express the ratio of 
some coefficients a, A in terms of T, once for (^l, and once for (/3r (see Appendix [B|1 : 



ORG 



2ie\/--l ^ 



2ieji 



1 



^1 

ARC 



■nil 



At. 



^-3>1 



< 1 



these inequalities become equalities only for T = 1. This proves that the symmetry 
between the regular and the singular part of a wave function which occurs at x = ±cxd 
is broken at a: = and that connection relations are not trivial (except for T = 1 and 
also the special case T = 0). 
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2.2.2. The S matrix The S matrix is related [TH [T5] to t and r and writes 

S-{: I). (15) 

Using the unitarity of the S matrix, it is useful to parametrize its elements in 
terms of the transmission coefficient T and a couple of two independent numbers 
€,e' = ±1. First, we get the parametrization of all coefficients ^l, ^r, which we 
give in Appendix [Bl Then, we can prove the representation 

^_ [T-l + iee WT-T ^ e'T + ieVT ^^ \ 

\ e'T + ieVT - T^ T-1 + ieeWT -T^) ^ > 

which is unitary, as required. We stress that this representation is not universal|16|. 
namely, it is peculiar to the Coulomb scattering problem as discussed here. 
We are now in a position to examine the connection problem. 



3. The connection problem 

The connection problem is to relate A, B to a,b either by finding matrix D in Eq. ([5]), 
or, equivalcntly, transfer matrix T in Eq. (|13p . or, equivalently, the S matrix in 
Eq. (|15p . Since ^ diverges as m ^ 0, it is not legitimate to use the continuity of 
(f and ^ at u = 0. Thus, the issue of the connection problem can not be handled in 
solving linear equations of the wave function, and one must address bilinear relations, 
related either to conservation laws or to certain constraints. In the following analysis, 
the behaviors of friiu), gri(u) and of their derivatives, for u ~ 0, are required: they 
are studied in Appendix [Cl 



3.1. Conservation laws and other constraints 



3.1.1. Continuity of p,, The simplest physical relation that provides a connection at 
x = is the continuity of the density of probability, Pr^iu) = [</'(«, r])\'^. With relations 
(|43al43bp . one gets 



-TTT] 



(17a) 



iBpe'^" = |6pe-"' ^ 

In Appendix |B1 we show that this relation actually simplifies as 

B = e'c-''%. (17b) 
where e' = ±1 (note that the case e' — —I implies a violation of the continuity of 



3.1.2. Current conservation The conservation of current j (x) = -3?(^iV,(x)^(x)j 
is equivalent to the unitarity of the S matrix which is already verified. Therefore, it 
does not help for the resolution of the connection problem. 



3.1.3. Orthonormality of scattering states Since the complete set of scattering wave 
functions is known, it is in principle possible to examine the consequence of generalized 
orthogonality relations. Let us write ip{x,E,a) = (pa{kx,jj:) with a = R, L (wave 
functions coming from +oo or — oo have degenerate energies), 

Jdx ■ip{x,Ei,ai)ii){x,E2,a2) = 5{ki - k2)Paia2 (18) 
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where P is an unitary 2x2 matrix in the (R,L) space. 

In Appendix |D1 using relations (|42al42bl42cl42dl42cl42t142g|42hp . ([T^ and 
(|12b|) . we calculate [TZj 



lim / ip{x, El, ai)ip(x, E2, a2)dx 

L—^oo 



xdiki -k2)+ \ + ee 



where c is a constant and Z a complex number given by, 

Z = V^(^i)^('72) + Vr(r?i)T(772) + iee'(Vi?fe)T(7n) - ^T{i^2)R{m)) ■ (20) 

Since fci,fc2 > here, we can drop 5{ki + fc2) in Eq. (fTO|) . which is irrelevant |18j. 
The established result in Eq. ([TO)) that P = I2 reflects the orthogonality of left and 
right moving states. Scattering states can be orthonormalized in the extended sense 
if and only if (/R^TX^ - ^ R{rii)T{r)2))Z = 0. This yields T{iu) = T{i^2) or 
T{r]i) G {0, 1}. The second condition is actually a particular case of the first one, 
since otherwise, one could find some energy E such that T{ri^) = 1 — T{ri^), which 
induces a non physical discontinuity; however, this argument will not be needed in 
the following. Having T independent of E is already a very strong result [TB]. Yet, in 
order to completely elucidate the connection problem, we will now address another 
constraint. 



3.1.4- Hermiticity of the Hamiltonian A successful issue for the connection problem 
is given by analyzing the hermiticity of Hamiltonian H. For Ei ^ E2, we consider two 
wave functions tpi : x ip{x, Ei) and 1P2 '■ x t-^ fp{x, E2) (degeneracy is not relevant 
here, and i?, L indices can be omitted). Since H is hermitian, the hermitian product 
of 1^1) with H\il)2) must be conjugate with the hermitian product of IV'2) with H\\pi). 
Explicitly, 



dxtp{x, El 



so that 



'^ix,E2)+^Mx,E2) 



-9V 



= dx 



ij{x,E2) 



dx ip{x,Ei)-^{x,E2) - -g^{x,Ei)ip{x,E2) = , 



'^{x,Ei)^{x,E2) + ^{x,Ei)iix,E2) 







(21) 



In Eq. (|2ip . we calculate the Cauchy principal value of the left term, which writes, in 
terms of dimcnsionless variables and function ip : 

-\ L 

>'/'(■", 772) 



lim — 

L — >QO 2 



ip{u,r)i) dip ^ dip , . . 
772 ou ou rji 



(22a) 
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Since —ip{u,i]i)^{u,7]2) + ^{u,rii)(p{u,T]2) is divergent at u = 0, one must use 
regularized integral around zero. Hence one should add the Cauchy principal value: 



lini — 



ipiu,7n) dip 



du 



— U,77l) 

ou ?]i 



(22b) 



and Eq. writes (|22ap + (|22bp =0. The contribution (|22ap is found to vanish when 
L ^ oo (detailed calculations, using relations (fT^ . ((i^ |42^ [i^gl [42h|) . are given m 
Appendix |E| so the net expression of Eq. pT|) is determined bv (|22bp which yields 



= Z 



3?(r(l+i772)-r(l+ir;i))VT(77i)T(772)} • 



Employing relations (|42al42bl42cl42d[) . we get the very same equation. Note that 



and h3{7]i,7]2) 



variable functions. Indeed, let us assume a linear combination. 



are independent two- 



Since 



jihi +72/12+73/13 = 
and 



(23) 



sinh(2:) 



are one-variable independent functions, if one keeps 772 



sinh(a:) 

constant and considers Eq. ((23|) as an equation of variable 771, one gets 71 = 0; if one 
keeps r]i constant and considers Eq. ([23|) as an equation of variable 772, one gets 72 = 0; 
thus, 73 = and the independence of the three functions is proved. Now Z, defined 
in (PO)) . can never vanish. Hence one gets 

R{ili)T{m) = ; T{7n)R{v2) = ; T{T^^)T{in) - . 



The first two equations imply T = 0, 1, and the last one simply implies T — Q. This 
eventually proves [T6] that, indeed, T(7^) = 0. 



3.2. Regularization by truncation 0} the potential 

Here we propose another approach, which gives the same result: the divergences are 
regularized by a truncation of the potential. 



3.2.1. Truncated half-potential In order to avoid the use of Coulomb wave functions 
for negative argument we calculate transmission and reflection amplitudes for a right 
half-barrier, defined for a; > 0, and then use reflection symmetry to calculate them 
for a mirror symmetric barrier, defined for x < 0. Then left and right barriers are 
combined using a composition formula for the S matrix, as suggested for instance in 
Ref. [IS]. 

The truncated right half-potential is, see Fig. [2l 

{0 for X < e ; 

A , (24) 
— tor X > e . 
X 

and the Schrodinger equation with Vs{x) alone writes — '^^^2^ +K(^)V'(^) = fc^-0(a;) . 

In order to avoid the \/\x\ singularity, the potential is assumed to be zero for 
< a; < e, but we have also performed our calculations with Ve{x<e) = -, with no 
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significant changes. The cutoff parameter e > is assumed small, and eventually the 
limit e — !■ is taken on the sum of left and right barriers, which corresponds to the 
complete Coulomb potential, since 




oJ I 

5 10 15 20 25 30 35 40 
X 



Figure 2. Right- half truncated potential II24I I and wave function in the two 
regions following Eq. ((26}. 

To calculate transmission and reflection amplitudes for the right barrier consider 
a plane wave approaching the potential 14 from — oo. It is partially reflected by 
the barrier at x = e, and the transmitted wave is a Coulomb wave tHrj, with 
Hri{u) = Fri{u) + iGn{u). Its asymptotic behavior is 



e 



i(u-e,,(M)) 



The scattering boundary conditions for the wave function are, see figure 

r ^ik(x-e) _^ -ik(x-e) for X < £ , 

I tHjjIykx) tor a; > e . 

We want to calculate reflection and transmission amplitudes r and t for this 
right-half truncated Coulomb barrier T4(a;). Matching a.t x — e yields 

l+r = tH,^(k£) , 1-r = -itH^{ke) , 

where H stands for dH/du, and thus 

i ^ 2 . ^ ^ H^{ke)+iH^{ke) 

if^ (fee) - iff,, (fee) ' Hriike) - iH,-,{ke) ' 

In the limit e ^ 0, this implies 

t ^ , r ^ -I . 

However, the limit e will not be taken here, but rather, at a later step. 

So far, we have considered transmission and reflection from the potential V^{x) 
where the incoming wave approaches the barrier from the left region. If the wave 
would have come from the right, be partially transmitted to the left and partially 
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Figure 3. Truncated potential for »; = 1, A = 1 and e = 1. 



reflected back to the right, the transmission amplitude would be the same, but the 
reflection would have a different phase. However, when we combine the symmetric 
image of (x) in order to account for the Coulomb problem as asserted in Eq. (pS)) , 
we employ the reflection amplitude r, as a result of the analysis developed in Ref. |19j . 
This procedure of combining the two barriers should be used before the limit e ^ 
is taken on Eqs. (|27p . The transmission amplitude through the combined barrier 

Ve{x) + Ve{-X) is 

J"^ = 

4 

= : : . (28) 

(1 - e2i'^^)[i?^(fce)2 - iJ^(fce)2)] - 2i(l + e2i'=-^)i?^(fce)77^(fce) 

This formula is exact and expresses the transmission amplitude for a symmetric 
combination of cutoff Coulomb barriers with a hole between —e and e. It uses 
Coulomb wave functions solely with positive argument. Inspecting the two terms 
of the denominator in Eq. (|28|) . the first term is found to vanish in the limit e —^ 0, 
and hence: 

i 

Te « : > . 

The upshot is that the transmission coefficient of combined left and right barriers, 
which comprise Coulomb barrier as e ^ 0, vanishes, that is T = VmiT^ ~ 0. 

3.2.2. A second form of truncated potential We also considered a truncated potential 
Ve, represented in Fig.[3]and defined as follow: e > and V|x| < e, Ve{x) = -, V|x| > e. 

The transmission can again be exactly calculated (the wave function 
corresponding to given {E, e) and its derivative are continuous; we use first order 
Taylor expansion for the Coulomb wave functions at connection points x = ±e). 

One finds, in Fig.[3]the curves of T!. versus e, for repulsive or attractive cases. We 
see that the transmission T!. ^ as e — *■ 0. This confirms our analytical result. We 
must precise that for some points of these figures, we used about 1000 digit precision 
calculation, provided by a formal calculation with integers. 



The one- dimensional Coulomb Problem 12 



0.0000 b 





Figure 4. Transmission versus e in the repulsive (plain line) or attractive 
(dashed line) case. 



4. Discrete spectrum : bound states 

We come now to the case of an attractive potential, and look for bound states of 
negative energies. As is shown below, analytical expressions can be obtained for the 
energies as well as for the wave functions |20|. 

4.I. Analytical solutions 

For e < 0, Eq. (|4]) is modified so that its right term writes ~(p(u) instead of (p{u). 
Note that u = kx holds but now k = \/— e, since, for an attractive potential, rj < 0. 
We will again consider separately u > and m < 0, and hence get the corresponding 
two equations: 

— ^(u) + 2-(, 
au u 



-j^{u) + 2-ip{u) ^ - ip{u) foru>0; (29a) 



, , (w) - 2-ip{u) = - ip{u) for u < . (29b) 
du^ u 

In order to solve Eq. (|29ap . we need to generalize equations (14.1.6), (14.1.14), 
(14.1.18), (14.1.19) and (14.1.20) of Ref. [10] (for L = 0). This is carried out in 
Appendix [Gl Generalization of (14.1.3) in Ref. [T^ is given below; relations (14.1.4), 
(14.1.5), (14.1.15), (14.1.17) remain valid by construction. Incidentally, the results 
of Appendix [G] can be regarded as a hyperbolic version of the original relations in 
Ref. [TU], since the solutions of Eq. (|29a[) now read: 

J^(u) = MC""il/(l + 77,2,2w) , A',,(u) EE 2ue-"/7(l + r?,2,2u) . 

In analogy with the case of free states, the functions and A'_^ are solutions of 
(j29bp (the connection problem at w = will be elucidated later on) . A useful identity, 
which will be needed, is 

J_„(U) = -J.r,{-U) . (30) 

^ . 2. Quantization 

For an arbitrary value of 77, the solutions Jr^{u) and Kjj{u) of Eq. (|29ap diverge as 
u ^ 00 and the solutions J-n{u) and K-^{u) of Eq. (|29bp diverge as u ^ —00. This 



The one- dimensional Coulomb Problem 



13 



is true for almost all values of 77, which therefore should be discarded as non physical, 
except for a set of quantized values ??„ (equivalently e„ or En) such that J^(u > 0) and 
J-ri{u < 0) are both square integrable, and for another set of values fjn (equivalently 
e„ or En) such that Kri{u > 0) and K-rj{u < 0) are both square integrable. The 
complete spectrum, which is described below, is composed of the union of set {En}, 
which is exactly Rydberg's spectrum, and set {En}-, the existence of which is indeed 
a surprise. 



4-2.1. The regular solutions Following the analysis of the hydrogen like atoms, it is 
verified that regular solutions J,^{u) and J-ri(u) decay exponentially as w — > ±00 only 
for a discrete set {?/„,Vn e N*} given by 

{qq')^m 

V = Vn = -n ^ E = En^-^^^—^. (31) 

The corresponding energies En form the Rydberg spectrum of hydrogen like atoms. In 
particular, the lowest energy is Ei = — 2(11^ yi^'i ~ —ZZ'Ei, where Ei is the Rydberg 
energy. 

The question whether the set 77,1 defined above can be used also for the singular 
solutions is answered negatively, although the demonstration is not immediate. While 
K-,f^(u) diverges as u —^ —00, K,f^(u) does not diverge as u — *■ 00. Therefore, one 
may consider a mixed solution AJri^ -\- BKrj^ for m > and aJ^ri„ for u < 0. However, 
as we shall see immediately below, Jt]„{0) ~ «^-?;„(0) = 0, while K-n{0^) = ^/C-nn- 
Hence the continuity of the density p at a; = implies here \B\ =0, which proves that 
a combination of regular and singular solutions is not an eigenstate. 

So far we have asserted the exponential decay of J±,,„ as u — ^ ±cx). The 
complete regular solutions Vn e N* can be constructed as Cn(u) = >^?;„(w) Vw > 
and Cn{u) = — M>^i)„(~") < 0, with £ C, (due to Eq. ([5D|) and the reflection 
symmetry between Eqs. (pgil) and (p9b)) ). Exphcitly (cf. Eq. (13.6.9) of Ref. [TU]), 

C„(.) = -^e-Hx;(2H)P J^'^^^^' (32) 
n \ 11 tor w < , 

where Ln{z) is the Laguerre polynomial of order n, and i^(^) ~ "^^dz^^ • ^^^^ 
shown below that /i = ±1. 

The orthogonality and normalization of the corresponding wave functions 
tjj{x,En) = Cnij^) ~ C'i('^^) inspected by carrying out integration on the 

positive semi axis K._|-. Thus, for the normalization we have, 

dx mx,En)f = / dx \Ukx)f du |C„(«)P " 



JO 



kj, '^"^ ^' |A|4 2|A| ' 



which, with \^\ = 1, requires a normalization factor equal to ^^^^ ; while for the 
orthogonality we find. 



dx tlj{x,En)i(j{x,En') = I '^"(^^)'^"'(^^) = yn^n' , 



due to orthogonality relations between Laguerre polynomials. 
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4-2.2. Anomalous solutions Quite remarkably, the anomalous solutions K^{u) and 
K_jj{u) both decay exponentially as u ^ ±00 only for a discrete set {77„,Vn G N} 
given by 

where £"„ is that of Eq. pip . The corresponding energies En form a separate spectrum 
interlacing the Rydberg one. From Eq. ([551) . one notes that £'„ — . j^i-.a ^'p; Vp e N*, 

SO that the minimum £^0 is lower than _Bi by a factor of 4. 

Note that, for rj 7^ f)n, K-,j(u) is diverging exponentially for u — > —00, while 
Krj{u) does not diverge for u 00. Therefore, one should examine the possibility of a 
continuous spectrum, by constructing a solution AKri{u) for u > and zero for u < 
for any such rj ^ ijn; however, one can calculate /ir,,(0+) = l/Cr, 7^ for all 77 < 0, so 
the continuity of the density p at x = implies A ~ 0. This possibility is eventually 
discarded. 

So far we have asserted the exponential decay of K±fi^ as u ±00. In order to 
construct the complete anomalous solutions, one needs to examine first the properties 
of K^fi^ (u) for u < and n G N. The imaginary part writes 

^{K-a (u)) ^Jf, {u) with 7„ {2n - 1)U/T+^ ; 

In 

while, for the real part, there is a relation analogous to ((30l) : 

/tt 

Kfi {-u) - i—Jfj (-m) = VnK-fj (u) \fu > (34) 

In 

where z/„ = 2^""''^/((2n + l){2n — 1)!!)^) : Kfj^ has even parity (whereas J,,„ has odd 
parity) if one omits rescaling factor z/„. 

The complete anomalous solutions Vn G N can then be defined as ^„(u) = Kfi^ (u) 
for u > and ^„(u) = vKf^^{~u) for w < 0, due to Eq. (j34p and the reflection 
symmetry between Eqs. ()29ap and (j29bp . It is not necessary to include the factor Vn 
here, since it is accounted for by the coefficient v. The latter will be shown below 
to be I' = ±1. In Appendix iHl we prove that the anomalous solutions are explicitly 
given by 

Uu) - (p„(H)Ko(H)+g„(H)Ki(|u|)) '"I X r II ' (35) 

(— 2)"Y7r I lor w < , 

where polynomials Pn{x) and ^^(a;) follow recurrence Eqs. ((i7a|) and ((47b)) . and K„ 
are the Bessel functions of the second kind. For instance, ^ qo ~ \, pi{x) = 3 — 4x, 
qi{x) = 1 — 4a;, P2{x) = 4x(4a; — 9) + 15 and P2{x) = Ax{4:X — 7) + 3 (more generally, 
these polynomials are proved to be real with integer coefficients in Appendix |H| . We 
are unaware of any occurrence of this family of polynomials, which are worth being 
studied further. 

As for determining the constant v, contrary to the regular case, ^„(0) 7^ 0. Hence, 
from the continuity of the density p, we deduce that 

16.(0-)! - |?„(0+)| 

in analogy with Eq. (|17ap . This implies v = ±1 (we are studying real solutions). 
Thus, the anomalous solution ^„ is even for z/ = 1 and odd for = — 1. 
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Similarly to the case of regular solutions, the orthogonality and normalization of 
the corresponding wave functions -(/'(a;, £'„) = ^„(^p) = ^„ ( J^Jj^j^ ) can be inspected by 
carrying out integration on the positive semi axis M+. Thus, for the normalization we 
have. 



JO 



2 ^ n u r m2 1 /(2n+l)/3„ , i^nT: 



dxmx,E„)\' = / dx\Ukx)\' ^ - / du\C„iu)\' - — 



kJo |A| V 22»+27r 2«+3 



The first coefficients /?„ can be easily computed, /3o = 3, /3i = 41, /32 = 1063. For large 
n, f3n ^ 5(2n + 1)!!. Since we proved v = ±1, one can deduce the exact normalization 
factor. 

Strikingly, the anomalous solutions are not orthogonal to each other. As a counter 
example, consider three hermitian products between anomalous states ^„ and with 
{n,p) = (0, 1); (0, 2) and (1, 2), on the semi-axis M+ : 



dxij{x,Eo)^{x,Ei)= I dxU\\\x)U^) 



Jo 



^ 3 _ 9(E(-8) - 3E(f) - 3K(-8) + K(|)) + 31n(729) 
IAPStt 64 

7^0.0210133 



J dxi;{x,Eo)iP{x,E2)^ J dxio{\\\x)^2C-^) 



^,35_ 175(E(-24) - 5E(||) - 4K(-24) + |K(§)) + 271n(5) . 
|A| ^ 487r 1728 

~ - 1^0.0319898 



Jlix^x,EiMx,E2)^ 

_ ^,i5_ _ 45(2705(3E(-f ) - 5E(f )) - 2877(5K(-f ) - 3K(f )) + 151n(729)) 
~ |Ap327r 256 ' 

~ T^O.0188906 

where K is the complete elliptic integral of the first kind and E is the complete elliptic 
integral of the second kind. It might be argued that these integrals were calculated on 
the semi- axis M+, while the hermitian product should be calculated on R and might 
vanish by symmetry cancellation (in case of odd parity, integrals on IR+ and on M_ 
have opposite sign). However, since we have already proved that all anomalous wave 
functions are either even or odd, then out of the three states (^Oj Cij C2), two have 
necessarily the same parity; thus, the corresponding scalar product is non zero, and 
these solutions are not orthogonal to each other. 

This is a surprising result which requires more insight into the properties of wave 
functions in quantum mechanics, which we will discuss briefly afterward. 



4-2.3. Orthogonality between regular and anomalous solutions Regular and 
anomalous solutions have different energies so they are expected to be mutually 
orthogonal as well (see also the discussion afterward). 
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Figure 5. Ci (red, full line), fa (purple, full line), §o (blue, dashed line) and ^2 
(magenta, dashed line). 



Performing the hermitian product on the semi-axis IR+ of -ip^x^En) = Cii(^^) 
with il>{x, Ep) = ( 2p+i ) yields a non zero result. For instance, 



^0 



dxil){x,Eo)il){x,Ei) = / dx£,Q{\\\x)Qi{ 



2 ^ 30F|A| 



similar expressions can be obtained for all n G N* and p g N, they can all be written as 
r/ ((7-y/7r|A|), with integers r and q depending on p and n. Thus, orthogonality between 
regular and anomalous wave functions can be assured only by symmetry cancellation 
of the right part of the hermitian product (on R-|_) with its left part (on K_). 

This leads to the following constraints: first, like the anomalous solutions, all 
regular solutions must have a definite parity. This is satisfied for /i = ±1. Second, all 
regular solutions must have the same parity, and all anomalous solutions must have 
the other parity. This means ^ = v is fixed. There remains a global choice of sign; 
either one chooses all regular solutions to be odd and all anomalous solutions to be 
even or vice versa. 

While we have no rigorous argument for either case, one notes that the choice 
^ = V = \ implies that C,'^ and ^„ are continuous. This seems to us the natural 
choice. Consequently, regular solutions Cn are odd and anomalous solutions ^„ are 
even. The first few solutions are shown in Fig. [5] With this choice, all solutions are 
continuous at m = 0, whereas the first and second derivative of ^„ are infinite at m = 
(this point is actually a ramification point). 



5. Discussion 



Despite its apparent simplicity, this one-dimensional problem leads to many interesting 
results, some of them unexpected. In the following, we will list our main results and 
discuss some of them. 



5.1. Zero transmission through the barrier 

The fact that T = for a repulsive infinite potential is in agreement with classical 
mechanics. On the contrary, for an attractive potential, it contradicts classical 
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mechanics. An example of perfect reflection from an attractive potential, called 
quantum reflection, is provided by the infinite square well potential: 

I 1 for bl < a , 
I for |a;| > a , 

where 2a is the width of the well. The Coulomb potential provides us with a 
new example of pure reflection. It differs from the infinite square well case by the 
width, which becomes narrower as one goes down in energy, and by the divergence of 
/ V{x)dx, which is logarithmic, while it is faster for the square well potential. Note 
that both the Coulomb potential and the infinite square well have an infinite number 
of bound states at negative energy. However, while the spectrum of the former is 
bounded from below, the spectrum of the latter is not. This is the only example of 
zero transmission and bounded spectrum that we know of. 

As a consequence of T = 0, singular unbound wave functions are eventually 
discarded, but the demonstration is much more involved than in the three-dimensional 
case of Eq. If one looks back at relations (|42al42bl42cl42dl42el42t142g|42hp . 

one finds that all B and b coefficients cancel: the logarithmic solution is completely 
suppressed, and therefore, the probability density is strictly zero at a; = 0. In the case 
of ipL, it is zero for a; > 0; in the case of i/j-r, it is zero for x < 0; the reflection process 
takes place entirely on the half-line. This suppression at a; = can be physically 
interpreted as a hard wall repulsion. It is also true for regular bound states, which 
probability density cancels at a; = 0. But it is not the case for anomalous bound 
states, which show, here again, a special behavior. 

5. 2. Representation of the S matrix 

In one-dimensional scattering problem with symmetric potential ^(a;) = V{—x), the 
S matrix is given by Eq. p5|) . Writing t = -s/Te'** and r = \J\ — Te'^'', the unitarity 
of S implies cos(6't — dr) = 0. Therefore, the most general expression of the S matrix 
can reduce to 



where < < Stt and e" = ±1 is an arbitrary sign. The form of Eq. (jl6p is unique to 
the Coulomb problem and reduces the number of free parameters, since there is only 
one continuous parameter T and two arbitrary signs e and e'. In particular, the phase 
dt is now given by 

tan(6'f) = , ^ . 

e'^T{\-T 

More precisely, (|16p relies on relations (|12ap . (jl2bp and on the reflection symmetry 
of the potential. For any symmetrical potential, one can choose a basis of solutions 
(/,(?) such that ()12ap holds; however, any generalization of relation ()17ap will fix the 
ratio A/a or B/b so that (|12bp will be changed. 

Note that, with T = 0, one simply gets S = —l2- 

5. 3. Non hermiticity of H 

The non orthogonality between anomalous bound states implies that H is not perfectly 
hermitian, because it is well established that the eigenstates of an hermitian operator 
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are orthogonal. This problem is raised by the same singularity than that, which is 
calculated in (|22bp . Indeed, the quantity A„p defined by 



2n+l 



2p- 



I ^^\x& 



2p+l, 



dx 



|A| 



/ 2n+l 



— lim 



d^p 

2p+l du 



2n + 1 



is not zero 



for instance Aqi 
etc. 



« A 



28 A 



12 



4 A 



23 
77r' 



23 14^ : 

But the situations are quite different. In the case of the unbound spectrum, 
eigenstates must be strictly orthogonal; otherwise, a quantum of a given energy E', 
coming from the frontiers of the universe and interacting with the system would not 
only create particles of the same energy, but of other energies, so E becomes blurred; 
but this blurring would spoil into the whole universe, which is impossible. So, we have 
discarded this possibility (proving therefore T = 0) of a break of hermiticity of H. 

On the other hand, a bound state of energy E may relax into a coherent state, 
thanks to interacting overlaps between non orthogonal eigenstates. Thus, it may be 
excited into a free state of different energy, with a certain probability, which we will 
examine; yet, this mechanism does not contradict any physical law, and is possible. 

Moreover, H is still an observable : its spectrum is reaZ, and canonical 
quantization theory is still valid, so a break of hermiticity restrictedly for E £ {E^ n S 
N} does not yield any contradiction of quantum mechanics, although its exceeds its 
standard axiomatic formulation. 



5.3.1. Coherent bound states Anomalous bound states are not orthogonal, so they 
are not stable: the spontaneous transition En Tp is allowed, without any interaction 
term in the Hamiltonian, which contradicts the standard properties of quantum 
mechanics. Therefore, a state of energy En is not stable. However, the transfer 
probability between two states of energies En and Ep is very small and decreases 
as \En — Ep\ is increased, so, anomalous states are almost stable, and their actual 
energy is only slightly blurred. In order to calculate stable states, one simply needs to 
diagonalize the (infinite) matrix M = {{S,m\^n))m.n- M is replaced by truncated matrix 
M'^^\ of size N X N corresponding to < m, rt < — 1, and we have diagonalized 
M'^^^ instead. By chance, the coefScients of M^^^ are rapidly converging when N is 
increased, so we can calculate numerically those of M . 

Let p(^) be the corresponding change of basis matrix, p(^) is indeed close to 
unity; we show, in Fig. [6] the rapid decrease of Pl'^fj versus N , in Fig. [7] the diagonal 
coefficient P^'^^ versus N, and, in Fig. [5] the convergence of P^^'' versus for some 

values of q (these coefficients are divided by P^"^^ for convenience). One verifies that 
the diagonal coefficient deviation from 1 remains very small, and, correspondingly, 
that other coefficients are of several orders smaller. 

The stable states that we have calculated are coherent states. Each coherent state 
can be labeled by the closest state of energy En and will be written When a state 
of energy En is created, it will relax to The delay of this relaxation is of the order 

, where AE'n is the uncertainty of En due to the instability process and can be 
explicitly calculated. 
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Figure 6. Pi'^ versus N. 



Figure 7. Af-[^' versus A' (it is normalized to 1). 




Figure 8. m{^' (red), m{^' (blue), M^^' (green) and m{^' (yellow) versus N 
(eoefBeient P^^' is divided by P^'''' to show tlie relative convergence). 



On the other hand, consider an excited state of energy E = ~Ep ; even if the state 
was initially created as 5„ with n ^ p, the probability of exciting state ^p, although 
small, is never zero. 



5.3.2. Orthogonality between regular and anomalous states Finally, we would like 
to insist on the orthogonality between regular and anomalous states. Otherwise, 
spontaneous relaxation between regular states. En Ep, might occur, through 
channel En — > Eq Ep, and the effective overlap between regular states would 
not be zero. 

If one adds, in the Hamiltonian, an interaction term between regular and 
anomalous term, allowing in-between transitions, the exact calculation of transfer 
probability would become more complicated, because of the relaxation process. 
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Eventually, in a real system, one should take into account the dynamical aspect 
of the problem, and consider, instead of a coherent state, an intermediate state, 
which would include the real dynamical relaxation process. Although it may seem 
complicated, this opens interecting directions. 



6. Conclusion 



Simple quantum mechanics can always bring new and surprising results. Indeed, we 
have found that the hermiticity of the Coulomb Hamiltonian may break exclusively 
for a closed family of bound states, which we therefore called anomalous states. These 
states are not stable, and one can only observe, instead, coherent states. We have 
also found a new case of quantum reflection, by solving the one dimension Coulomb 
problem. 
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Appendix 

A. Asymptotic behavior of F^i and Gn when u — > —oo 

Here we analyze the asymptotic behavior of F,^{u) and Gri{u) for u — > — oo. Our 
results are different from those in Eqs. (6c) and (6d) in Ref. [I], (see |13)). 
Let us first demonstrate (|llal) . First note that 



ie**M(l + i?7,2,-2iu) = tc-"M(l - iT],2,2iu) , (36) 
but, since it is real, one can omit the conjugation. For u > 0, writing u = \u\ and 



using (|9a|) . one gets 

|M|e~*l"lAf(l-i77,2,2i|w|)|jp^e'?K^sin(|u| -e^(u)) . 
For u < 0, writing u ~ — |u| and using ((36)) . one gets 

|u|e*'"lM(l + i7;,2,-2i|u|),'^ c'^K^jSindul - e^(u)) ; 

V 

— — sin(ii+0T7(u)) 

if you make 77 —77 in the last relation, and multiply by -1, you get 

-\u\ e^l'^lM(l - 2, -2i|w|)|^^p-j^e~^K^ sin(ii - B^(w)) , 
which is exactly the expected relation 

^'7(-l"l)i„^i:T;oC""''sin(^.-e,(j.)) . 

We only write here the leading order of (jllap . you must be very careful of all sign 
compensations for the next orders. Eventually, if one makes again 77 —* —77 in the last 
relation, one gets directly 

F„,,(-|u|) -~ e-''sin(7. + e^(7.)) , 
which is the behavior of jr]{u) for u ^ —oo. 
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The demonstration is very similar, for (jllb|) . First note that 



te"*C/(l + ir;,2,-2iii) = ic-"J7(l - i?7,2,2iM) ; (37) 

here, conjugation can not be omitted. For u > 0, writing u — \u\, using (j9bp and 
keeping only the real part, one gets 



3? 



f |.| e-*l"lt/(l - i,, 2, 2i|u|)) , ~ j; ^o<\-\-^niu)) 



l"H+°o2ry3?(r(-i?/)) 

For u <0, writing u = — \u\, using (|37|1 and still keeping only the real part, one gets 
^ (l"l e-'"l Uil + i,, 2, -2i|.|)) |.~oo^^;^4^ • 

— cos('(iH-0,, (u)) 

if you make — r/ in the last relation, and multiply by -1, you get 

e^ cos{u — Qriiu)) 
l'^P+^2r/3fi(r(-iry)) ' 



3? (^-|M|e*l"li7(l-i7?,2,-2i 
which is exactly 



Eventually, if one makes again r/ —t] in the last relation, one gets directly 
G_,(-H)|~^c--''cos(^ + e,(u)) , 

which is the behavior of gr]{u) for u — oo. 
B. Expression of t as a function of T 

First, you get simple relations between {ta,ra) and {Aa, Ba,aa,ba) (a =R,L) : 

(38a) 
(38b) 
(38c) 
(38d) 
(38e) 
(38f) 
(38g) 
(38h) 

The unitarity of S writes 

|r|2 + ^ 1 ; (39a) 
tr + rt = . (39b) 





= i^L ; 






= iL ; 






= ie-'^"! 




6l 


= e'^''(l 


+ 7-l) ; 




= -i(l 


- m) ; 


Br 




7 


Or 


= — ie~ 


""tR ; 


bn 







From (j39bp one deduce 
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where e = ±1. From relations (j38bl38dp . one gets 

By use of relations (|39a|39bp , this writes 



5l t 
but Eq. (|17ap implies the existence of 6* e M such that 

so, using back relation p4)) . we get 



i^e"-i,^i,-l. 

We carefully multiply this equation by its conjugate and find 



which implies = or tt. We will write e'^ = e' then 



We carefully multiply this equation by its conjugate and find 
but |tp = 5R(t)2 + so we get 



By use of (0, we have t = ^{t) + = e'T + ie'VT - T^. We eventually shall 

prove that e" = e. We put the last expression of t into (1 + r)/t and get 



1 + ^ 1 + iety^i - 1 (1 - £e" + re(e" + ie'y^i - l))(e'T - ie'VT - T2) 



t 



By taking the modulus of this expression, you would find indeed that e = e". However, 
we already know that it is real (because = or tt), so you have the result straight. 
Now, if you use back the different relations, you can get the final expression of T : 



t = e'T + ieA/T(l - T) (41) 
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where e' = ±1 is independent of e. By use of relations (|38el38f138g|38iil38al38 bl38cl38d|) . 
TTj) and (|17ap . after some calculations, one gets 



e'T + ie^T{l-T) ; 
e-''"(eeVr(l - T) + i(2 - T)) 
e"''(T + ieeVr(l-r)) ; 



= e-'^''(eVT(l - T) - ie'T) 
e^^ie'T + ieVTil-T)) ; 



Ar = - eeWTil-T)- i(2 - T) 
Br, = T + ieeVni - r) ; 
Or 



and 



r = T - 1 + iee'y/T{l-T) . 

Using these relations, one verifies all relations (|14l39al39bp and (jl7bp . 

An important collateral result from this demonstration is indeed that 



&Le~ 



from relations (j38fl38bl38hl38dp . one gets 

Bn ~ Bl 
which proves, by linearity, relation (|17bp . 



(42a) 
(42b) 
(42c) 
(42d) 
(42e) 
(42f) 
(42g) 
(42h) 

(42i) 



C. Mclaurin expansions 

Here we study the behavior of basic solutions f.q{u), gri{u) and their derivatives when 
u — > 0. Let us consider first the expansion of F^j and G,, for u ^ 0^ , which are given 

hym ■■ 

Frj{u) ~ C-T|r(l+i77)|(u + ?/t2) 

G,,iu) ^ !^2rjiu + r,u^) ( log(2w) - 1 + p{ri) + 27e) + (1 - ^-^p^u^) 
au 

dG„ , , 1 
-(^) - 7^ 

■n 



-^(w) ~ — {2?7 [(1 + 2r;u)(log(2u) +p(r;) + 27e) - - (1 + 



du 
d^F 



dv? 
d^G 



277 ( log(2u) + p{7]) + 27E ) + ?? + - 



- (1 + 677' 
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with p{i]) = 3? ( ^r(i+i!,) ) = Pi^v) ^i^d 7e is Euler's constant. Thus, one gets, at first 
order, for the complete solution </?, 

v[n,Ti) -^_b = 7^; (43b) 

^{u,in)-~^^AC.^ + 2B'n^{\og{2u)+p{'n) + 27E) ; (43c) 



dip,^ ^ ^ 1 

-n 



^{u,Tj)^_aC^^ - 2brj—i\og{^2u)+pirj) + 2^e) ■ (43d) 



D. Orthonormality relations 



The purpose of this section is to calculate the limit, when L oo oi 

J^j^'ip{x, Ei,ai)'ip{x, E2,Oi2)dx. Consider a given L, this integral with all functions 



replaced by their asymptote (|12a|) or (jl2bp becomes: 



dx cos 



cos 



\x 



dx cos 



cos 



2 Vim 

r}i-ri2 

Xx 

2 Vll2 

vi+'n2 
Xx 

V1+V2 

Xx 
2 Vim 

V1-V2 

Xx 
2 mm 

VI -'V2 

Xx 



xX 



0r,.( 



2r]i 

0m(i^)-e..(7F) 



xX 

e,.(|^) 

xX 



xX 



0>u(;^) 

2r/i 
xX 



2772 

Ov2i^) 



xX . 



a; A 
2^ 



4 + 



A- 



sm 



sm 



-) V1V2 

Aa; 

> VlV2 

' V1+V2 

Xx 

> VlV2 

'Vi—V2 



+ Qm(:r-) + e,J— 



27? 



where we use 



A + 

aiQ2 
aiQ2 



b+ 



CKl a2 



^Q2 ^ 



■7 7r(??i-^?2) I U 



a, 



,'r(')l-f;2) 



^Qia2 (^ai-^a2 -^01^02) 7 
-7r(?)i+r;2) ^ 

-7r(?)i+T;2) ^ 

-7r(»;i-r;2)-j 

-7r(r;i-T;2)\ 
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The difTcrcncc with the exact hmit is finite and contributes to constant c in 
formula p9p . Now. these integrations are easily performed when one notes that 
all Sr,(u) functions can be treated as constant. Indeed, let us consider a simpler 
integral J^^ cos{su + ln{u))du, where we will omit the problem at u = 0, and 
6{L) = -J sin(su + ln(w)) — J^^ cos(su + \ii(u))du is the difference of the approximate 
integral with the exact one. Then, S'{L) = ^'"(^^+^'"(^)) not only tends to zero when 

L ^ oo, but has a finite integral 5'{u)du. This proves that all such approximations 
are valid and simply contribute to constant c. 

The .T = boundary only contributes to constant c (you may need to replace 
X = Q with another boundary, in order to avoid any divergence, but this replacement 
simply gives another contribution to constant c) so we may skip it and eventually get 



Vim 

m - m 
mm . 



■ sm 



+ 



m + m 

mm 
m + m 

mm 
m - m 

mm 

m - m 

mm 
m + m 

mm 
m + m 

mm 
m - m 



■ cos 



■ COS 



■ COS 




Now, both limits of 



sin(sL) 



171 -i;2 

and 22£(£M 



when L — > 00 are equal to t:5{s) (with 



differential ds). The ln(u) correction has no influence (see Appendix [E| . Then we 
i-) = S{\{ki — k2)) = ^5{ki — fc2), so we eventually get factor We 



write 5( — 

have forgotten the exact differential ^ in one dimension, and we will include a last 

factor 2 which accounts for the equality between the limits of and /_^- Altogether, 
we get formula (fT9|) , with the following coefficients of matrix P : 



Pr. 



An An 



BaBa' + OaOa' Q^'^'^ + hobo 



and, with relations (|42el42t142g|42hl42al42bl42cl42dp . we eventually get 



P 



thus (flSl) is verified. 
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E. Hermiticity relations at infinity 



The calculation of (|22ap is similar to the previous orthonormality calculations, 
although simpler. Here a = R, L for the choice of ipa and we use the notations 
of Appendix |DJ One gets 



m + m 
vi - m 

27?1??2 

T?i - m 

2r]ir]2 

2771772 
vi + m 

2f7i'72 

^ vi - m 
2771772 



XL 



) mm 
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sm 
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2771772 
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;> 1)11)2 
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;> ')i')2 

'^vi-m 
XL 

) 1)11)2 
" 1)1 -1)2 

AL 

1 1)11)2 

' 1)1+1)2 

AL 

-> 1)11)2 
■^»)1+1)2 

AL 

2 1)11)2 
■^1)1-1)2 



*»'^'- 




. _|_ 






1 ^aa 


2r/i 


r. ,LX 
27?2 


1 / 1 1 


LA 




1 B^a 






) 


2771 












''''^4771^ 


,LX 

+ 01)2^^ 

27?2 


1 ^aa 



One striking thing is that the coefficients — ± — are very different from the 

00 1)1 1)2 

previous case. In order to match with the 5 limit, one must divide by T so there 

is a global supplementary factor ~ which, when multiphed by (5(^ ± will 
always give zero. 

Another important difference is that we have made no approximation in this case. 
It is worth studying the last limit more carefully, than we did before. Using again a 
simpler case, we want to prove that limi^oo 7 sin(sL — sln(L) — ks^ + /3) is 7r(5(s) (k 
and (3 are just constants here). The important thing is that L = L — ln(L) 00 and 
can be used as a parameter, so the result is proved, and the limit of (j22ap is strictly 
zero. 



F. Digression: To WKB or not to WKB? 

In a nuclear fission process, a light nucleus of mass tti and charge q = Zqe > (c.g 
an alpha particle with Z = 2) is trapped in a metastable state at energy E due to 
a potential "pocket" V{r) = V^Ar(r) + Vc{r) of an heavy nucleus of charge q' = Z'q^ 
(here r is the distance between the centers of mass of the two nuclei). The potential 
is the sum of a strong short-range attractive nuclear potential Vi\f{r) and a repulsive 
long range Coulomb potential Vdr) ~ The focus of interest is on the escape 
probability P from the metastable state. In a crude approximation, V{r) is replaced 
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by a deep potential well of range R and depth — Vq and a Coulomb tail for r > R, see 
figure [5) 




Figure 9. Approximated potential V{r) = V]\f{r) + Vo(r) designed to calculate 
the fission probability within the WKB approximation. 

The escape probability is then calculated in the WKB approximation, integrating 
the local momentum K(r) = ^ 2™ ^y^ (7.) _ between the turning points R and Rc 
(such that Vc{Rc) = E). 



Ar 



A. 



~2A 



/•flc /--Rc rR 

A = / K{r)dr = / K{r)dr~ I Hi{r)dr 
Jr Jo Jo 

When R Rc the result is written as, 

27vmqq' 32mqq'R 

p = e — e r^'^ = PgTr, (44) 

where v is the relative velocity and Pq is the Gamow factor, which contains the energy 
dependence of the escape probability. Relation ([27|) , with e — R, gives the exact escape 
amplitude ~ \t\'^ (for the special case Vb = but that can be easily modified). It also 
shows that the WKB expression (|44p cannot be used as i? ^ because it yields a 
finite escape probability while the exact result (within the naive model of figure [5]) 
gives zero escape probability. The reason is that the conditions for the use of the 
WKB approximation arc not met, strictly speaking. 



G. Generalization of recurrence equations 

We study the changes of relations (14.1) in Ref. [TD| for the bound states (e < 0), in 
the case L = 0. Note first that (14.1.1) is also changed, it writes now as (|29ap . 
Relation (14.1.6) writes now (we omit the L = exponent) 

Ai = l; A2=ir, (fc + l)(fc + 2)Afc+2 = 2?7^fc+i + Afc . 

Relation (14.1.14) writes now (with our notations) 

L^{u) = 2ryX,(ii)(log(2u) - 1 + + 27e) + e.,{u) 

with (14.1.17) (relation (14.1.15) is useless here) 

oo 
fe=0 
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and relations (14.1.18) to (14.1.20) now become 

ao = 1 ; ai = -1 ; (fc + l)(fc + 2)ak+2 = '^vo-k+i + Ofe - 2ij{2k + 3)Ak+2 ■ 

Eventually, note that new relation (14.1.14) also holds for u < as soon as you replace 
log(2w) by log(-2u). 

H. Identities between confluent hypergeometric functions an modified Bessel ones 

We found useful identities between confluent hypergeometric functions M(i ±n, 2, 2t) 
or J7(i ± n, 2, 2t) and modified Bessel functions I„(t) or K„(t), for all n E N. 

These identities appear to generalize some identity established only for n = or 
n = 1 ; indeed, from relations (13.6.3) and (13.6.21) of Ref. [10|, one shows 

e-*A/(i,2,2t) = Io(t)-Ii(0 ; (45a) 

e'*C/(i, 2, 2t) = ^(Ko(t) + Ki(t)) ; (45b) 

Thus, it seems possible to generalize this relations and look for solutions of Eqs. (|29ap 
and (|29bp in the form 

fn{t) = t(p„{t)Mt) - (?n(i)li(0) (46a) 

or 

gn{t) = t{p„{t)K„{t) + <7„(t)Ki(0) (46b) 

(we took advantage of further relations between the polynomials {pm<ln) defined in 
Eq. (|46ap and those defined in Eq. (j46b|l in order to save notations.) 

Although it works well, it proved more efficient to find directly the recurrence 
relations which define p„ and (j„. Using relation (13.4.11) of Ref. [TOj for AI{^—n, 2,2t), 
(13.4.10) for A/(i+n,2,2t), (13.4.26) for U{\~n,2,2t) or (13.4.23) for C/(i + n, 2, 2<), 
and making the derivative of Eqs. (|46ap and (j46bp using Ig = Ii, Ii(i) = Io(i) — Ii(t)/i, 
Kq = — Ki and K']^(t) = — Ko(t) — Ki(t)/t, and fixing Po = go = li one finds, up to 
some normalisation factors, 

Pn+i{x) = (2n + 'i)pn{x) + 2a;(p^(x) - _p„(x) - (7„(x)) (47a) 
(7„+i(a;) = (2n + 1)(?„(.t) + 2x(q'^(x) -_p„(a;) - g„(a;)) (47b) 
These definitions have one main advantage : these polynomials are real and have 
integer coefficients ; let us write Pn{x) = X]r=o ^"^^ <ln{x) = X]"=o^?^'' '^^ 
= (2n + 1)!!, b^ = {2n - I)!!, al = bl = (-4)". 
Eventually, let us fix the normalisation problem (note the symmetry between 
M{\ - n, 2, 2t) and M(| + n, 2, 2t) or between [/(i - n, 2, 2t) and f/(| + n, 2, 2t) and 



that (-1)!! = 1): V?i € N, 

G-'M{\ - n,2,2t) - ^--l-^(p„(t)/o(t) _ g„(t)/i(0) ; (48a) 

1 C— 1")" / \ 

e-*C/(- - n,2,2t) = A^^(p„(t)A-o(t) +g„(t)i^i(t)) ; (48b) 

e-*M(^ + n,2,2t) = — -1— ^(p„(-t)/o(i) + g„(-t)/i(t)) ; (48c) 



3 2" / \ 

e-C/(-+.,2,20 - (2n+l)!!(2n-l)!!V^ (-^"(-^)^"(^)+ 
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